Topological order in the insulating Josephson junction array. 
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We propose a Josephson junction array which can be tuned into an unconventional insulating 
state by varying external magnetic field. This insulating state retains a gap to half vortices; as a 
consequence, such array with non-trivial global geometry exhibits a ground state degeneracy. This 
degeneracy is protected from the effects of external noise. We compute the gaps separating higher 
energy states from the degenerate ground state and we discuss experiments probing the unusual 
properties of this insulator. 



A device implementing an ideal quantum computer 
would be a very interesting object from a physics point 
of view: it is a system with exponentially many degen- 
erate states (TV ~2 X , K ~ 10 4 ) and extremely low de- 
coherence rate 1 . The latter implies that the coupling 
of these states to the noise field produced by environ- 
ment should be very small, more precisely all states 
are coupled to the external noise in exactly the same 



way: (n 



O 



m) — OoS mn + o(e L ) where O is an op- 
erator of the physical noise and L is a parameter such 
as the system size that can be made as large as de- 
sired. This condition is naturally satisfied 2 if degenerate 
states are distinguished only by a non-local topological 
order parameter. 3 ' 4 . In the previous work 5 we have pro- 
posed a superconducting Josephson junction array that, 
in addition to superconducting long range order, acquires 
a topological order parameter for the arrays with non- 
trivial geometry. The goal of this paper is to identify the 
Josephson junction array with the ground state charac- 
terized only by the topological order parameter and no 
other long range order of any other kind (in local vari- 
ables) and which has a gap to all other (non-topological) 
excitations. Because such array has neither supercon- 
ducting phase stiffness nor gapless excitations it can not 
carry an electric current; thus in the following we shall 
call it topological insulator. While the arrays with simi- 
lar global properties were proposed and studied in a re- 
cent work 6 , the arrays discussed here have a number of 
important practical advantages which makes them more 
feasible to built and measure in laboratory. 

In physical terms, the topological superconductor ap- 
pearing in the array 5 is a superfiuid of 4e composite ob- 
jects. The topological degeneracy of the ground state 
arises because 2e excitations have a gap. Indeed, in 
such system with the geometry of an annulus, one ex- 
tra Cooper pair injected at the inner boundary can never 
escape it; on the other hand, it is clear that two states 
differing by the parity of the number of Cooper pairs at 
the boundary are practically indistinguishable by a local 
measurement. 

Generally, increasing the charging energy in a Joseph- 
son junction array makes it an insulator. This transition 



is due to an increase of phase fluctuations in the origi- 
nal array and the resulting appearance of free vortices 
that form a superfiuid of their own. The new situa- 
tion arises in topological superconductor because it al- 
lows half-vortices. Two scenarios are now possible. The 
"conventional" scenario would involve condensation of 
half-vortices since they are conjugate to 4e charges. In 
this case we get an insulator with elementary excitations 
carrying charge 4e. An alternative is condensation of 
full vortices (pairs of half vortices) with a finite gap to 
half vortices. In this case the elementary excitations are 
charge 2e objects. Similar fractionalization was discussed 
in the context of high T c superconductors in ,s and in the 
context of spin or quantum dimer systems in 9-13 . Such 
insulator acquires interesting topological properties on a 
lattice with holes because each hole leads to a new binary 
degree of freedom which describes presence or absence of 
half vortex. The energies of these states are equal up to 
corrections which vanish exponentially with the size of 
the holes. These states cannot be distinguished by lo- 
cal measurements and have all properties expected for a 
topological insulator. They can be measured, however, 
if the system is adiabatically brought into the supercon- 
ductive state by changing some controlling parameter. In 
this paper we propose a modification of the array 5 that 
provides such control parameter and, at the same time, 
allows us to solve the model and compute the properties 
of the topological insulator. 

The main physical idea of the array is to give a system 
where kinetic energy of half vortex, thv is parametrically 
smaller than the kinetic energy of the full vortex, tf v and 
where their potential energies, W, satisfy th v < W < tf v . 
The array is shown in Fig. 1, it contains rhombi with 
junctions characterized by Josephson and charging ener- 
gies Ej > Ec and weak junctions with ej ec <C Ec- 
Each rhombus encloses half of a flux quantum leading to 
an exact degeneracy between the two states of opposite 
chirality of the circulating current 5,14 . This degeneracy is 
a consequence of the symmetry operation which combines 
the reflection about the long diagonal of the rhombus and 
a gauge transformation needed to compensate the change 
of the flux $o/2 — > —$0/2. This gauge transformation 
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changes the phase difference aiong the diagonal by ir. 
This Z 2 symmetry implies the conservation of the par- 
ity of the number of pairs at each site of the hexagonal 
lattice and is the origin of the Cooper pair binding. We 
assume that each elementary hexagon contains exactly k 
weak junctions: in case each link contains one weak junc- 
tion k = 6, but generally it can take any value k > 1. 
As will be shown below, the important condition is the 
number of weak junctions that one needs to cross in the 
elementary loop. Qualitatively, a value k > 1 ensures 
that it costs a little to put vortex in any hexagon. 

For the general arguments that follow below the actual 
construction of the weak links is not important, however, 
practically it is difficult to vary the ratio of the capaci- 
tance to the Josephson energy so weaker Josephson con- 
tact usually implies larger Coulomb energy. This can be 
avoided if weak contact is made from Josephson junction 
loop frustrated by magnetic field. The charging energy of 
this system is half the charging energy of the individual 
junction while the effective Josephson junction strength 
is ej = 2tt^Eq where E Q is the Josephson energy of each 

contact and 5$ = $ — % is the difference of the flux from 
half flux quanta. This construction also allows to control 
the system by varying the magnetic field. 

Assume that ej sets the lowest energy scale in this 
problem (the exact condition will be discussed below). 
The state of the array is controlled by discrete variables 
u a b = 0,1 which describe the chiral state of each rhombus 
and by continuous phases cf> a b that specify the state of 
each weak link (here and below a, b denote the sites of 
hexagonal lattice). If Josephson coupling ej = different 
islands are completely decoupled and potential energy 
does not depend on discrete variables u a t. For small ej 
we can evaluate its effect in the perturbation theory: 

V(u) = -Wcos(irY,u ab ),W=-(-^\ ej (1) 

hex ^ ^ 

This potential energy lowers the energy of classical con- 
figurations of u a b that satisfy the constraint ^2 hex u = 
[2] but does not prohibit the ones with Y^hex u = •"■ P]- 
Consider now the dynamics of discrete variables. Gen- 
erally, two types of tunneling processes are possible. In 
the first type the phase changes by ir across each of the 
three rhombi that have a common site. This is the same 
process that gives the leading contribution to the dynam- 
ics of the superconducting array 5 , its amplitude is given 
(in the quasiclassical approximation) by 

r «££ /4 ££ /4 cxp(-3S ), S = (2) 

In the second type of process the phase changes across 
one rhombus and across one weak junction. Because 
the potential energy of the weak junction is assumed 
to be very small the main effect of the weak junction 
is to change the kinetic energy. The total kinetic en- 
ergy for this process is the sum of the terms due to the 
phase across the rhombi and across the weak link. As- 




FIG. 1: Schematics of the array. The main figure: global 
structure of the array. Discrete variables controlling the low 
energy properties are defined on the links of the hexagonal 
lattice. Generally, the lattice might have K big holes, here 
we show K — 1 example. Zoom in: Each inner bond of the 
lattice contains a rhombus made out of four Josephson junc- 
tions; some bonds also contain an effective weaker link made 
of two Josephson junctions so that each hexagon of the lattice 
contains k = 3 such links. The flux through each rhombus is 
half-flux quanta, 3>o/2, the flux through a loop constituting 
a weak link is close to half flux quanta <E> = $o/2 + 5$. The 
boundary of the lattice contains rhombi and weak links so 
that each boundary plaquette has the same number, k, of 
weak links as the bulk hexagon. 



suming that these phase variations are equal and oppo- 
site in sign, the former is about Eq (j) 2 while the latter 
e^ 1 ^ 2 , so the effective charging energy of this process is 

E c = (Eg 1 + e^ 1 ) . For e c <C E c this charging energy 
is small and such process is suppressed. Thus, in these 
conditions the dominating process is the simultaneous 
flip of three rhombi as in the superconducting case. In 
the following we restrict ourselves to this case. Further, 
we shall assume that r > IV so that in the leading or- 
der one can neglect the potential energy compared to the 
kinetic energy corresponding to the flip of three rhombi. 
As W is increased by turning on ej the continuous phase 
4> a b orders and the transition into the superconducting 
state happens at ej ~ ec- At larger ej, W becomes ej 
and with a further increase of ej, for ej 3> r vortices 
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completely disappear from the low energy spectrum and 
the array becomes equivalent to the one studied in 5 . 

The low energy states are the ones that minimize the 
kinetic energy corresponding to simultaneous flip pro- 
cesses: 

H T = -rJ2U T « b ( 3 ) 

a 6(a) 

Here b(a) denote the nearest neighbors of site a, r x b is the 
operator that flips discrete variables u a b and r is given by 
(2). The states minimizing this energy satisfy the gauge 
invariance condition 

n^i*)=i*> ( 4 ) 

6(a) 

The Hilbert space of states that satisfy the condition (4) 
is still huge. If all weaker terms in the Hamiltonian are 
neglected all states that satisfy (4) are degenerate. These 
states can be visualized in terms of half vortices posi- 
tioned on the sites of the dual lattice, Indeed, a 
convenient way to describe different states that satisfy 
(4) is to note that operator Ylbta) T ab d° cs not change 
the value of X)j(i) u P]- Thus, one can fix the values of 
u ij = v i on au plaquettes, i and impose the con- 
straint (4). In physical terms the binary values Vi = 0, 1 
describe the positions of half-vortices on dual lattice. 
This degeneracy between different states is lifted when 
the subdominant terms are taken into account. The main 
contribution to the potential energy of these half- vortices 
comes from (1), it is simply proportional to their number. 
The dynamics of these vortices is due to the processes in 
which only one rhombus changes its state and the cor- 
responding flip of the phase accross the weak junction. 
The amplitude of this process is 

r » Ey^E^expi-So), So = 1.61 J 

V Ec 

The effective Hamiltonian of half- vortices is 

(ij) » 

where operators Ui act in the usual way on the states 
with/without half- vortices at plaquette i and the first 
sum runs over adjacent plaquettes (ij). This Hamiltnian 
describes an Ising model in a transverse field. For small 
W/r < A c <~ 1 its ground state is "disordered": (a z ) = 
but (<j x ) ^ while for W/r > A c it is "ordered": (a z ) ^ 
0, (a x ) — 0. The critical value of transverse field is known 
from extensive numerical simulations 15 : A c « 4.6±0.3 for 
triangular lattice. The " disordered" state corresponds to 
the liquid of half- vortices, while in the "ordered" state 
the density of free half-vortices vanishes, i.e. the ground 
state contains even number of half-vortices so the total 
vorticity of the system is zero. To prove this we start from 
the state | j) which is the ground state at r/W — and 



consider the effect of r X^(jj) a j a f m perturbation theory. 
Higher energy states are separated from the ground state 
by the gap W so each order is finite. Further in each order 
operator a x af create two more half-vortices proving that 
the total number of half-vortices remains even in each 
order. The states with odd number of half-vortices have 
a gap A(r/W) which remains non-zero for W/r > A c . 

In terms of the original discrete variable defined on the 
rhombi the Hamiltonian (5) becomes 

^ = -^E^-^£IK- ( 6 ) 

(ij) 1 }(*) 

where operators act on the state of each rhombus. This 
Hamiltonian commutes with the constraint (4) and is in 
fact the simplest Hamiltonian of the lattice gauge 
theory. The disordered regime corresponds to a confined 
phase of this Zi gauge theory, leading to elementary 4e 
charge excitations and the ordered regime to the decon- 
fincd phase. 

Consider now the system with non-trivial topology, e.g. 
a hole. In this case the set of variables Vi is not sufficient 
to determine uniquely the state of the system, one has 
to supplement it by the variable vq — Y] L Uij where sum 
is taken over a closed contour L that goes around the 
hole. Physically, it describes the presence/absence of the 
half-vortex in the hole. The effective Hamiltonian of this 
additional variable has only kinetic part because presence 
or absence of half vortex in a hole which has I weak links 

in its perimeter gives potential energy Wo = cej (f^) 
which is exponentially small for I 3> 1. The kinetic part 
is similar to other variables: Ho = — r^2 ieI erf crfi , it de- 
scribes a process in which half- vortex jumps from the hole 
into the inner boundary, /, of the system. In the state 
with ((T z ) 7^ this process increases the energy of the 
system by W(^) (W(0) = W and W(\ c ) = 0). In the 
state with (<j x ) ^ it costs nothing. Thus, the process 
in which half- vortex jumps from the hole into the system 
and another half-vortex exits into the outside region ap- 
pears in the second order of the perturbation theory, the 
amplitude of this process is t v — r 2 ^2 ieI je0 gij where 
sum is performed over all sites of the inner (J) and outer 
(O) boundaries and = (af H l Eo &j^ has a physical 
meaning of the half vortex tunneling amplitude from in- 
ner to outer boundaries. At small r/W, we can estimate 
gij using the perturbation expansion in r/W: the lead- 
ing contribution appears in \i — j|-th order of the per- 
turbation theory, thus oc (r/W)^^ . Thus for small 
r/W the tunneling amplitude of the half vortex is expo- 
nentially small in the distance, L, from the outer to the 
inner boundary; we expect that it remains exponentially 
small for all r/W < X c . For r/W > A c this amplitude is 
of the order of r 2 /W and therefore is significant. 

In a different language, in the system with a hole we 
can construct a topological invariant V = J7 (con- 
tour 7 is shown in Fig. 1) which can take values ±1 . 
The same arguments as used for the superconducting ar- 
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ray show that any dynamics consistent with constraint 
(4) preserves V. Thus, formally the properties of the 
topological insulator are very similar to the properties 



of the topological superconductor discussed in 5 , if one 
replaces the words Cooper pair by half-vortex and vice 
versa. We summarize this duality in the following table 



Topological Superconductor Topological Insulator 

Ground state Condensate of 4e charges Condensate of 2ir phase vortices 

Fluxons Capful, charge 2e Gapful, tt phase vortices 

Pseudocharges Half fluxes with energy e ~ Ej log(L) Charge 2e with e = 2r. 17 

Ground state degeneracy Charge on the inner boundary mod 4e Number of ir vortices inside the hole mod 2 

Ground state splitting (^^) L (r/W) L 



Note that at small r/W — ► the ground state of 
the Hamiltonian (6) satisfies the condition (4) and mini- 
mizes the second term in (6), i.e. satisfies the condition 
Ilj(i) T ij\^) — 1^); it can be explicitly written as |0) = 

Il< U 1 + Uj(i) Tij)Xlkl I -»>*!■ This state is a linear su- 
perposition ^(|+) + | - )) °f the degenerates states with 
V = 1 and V = —1; it coincides with the ground state of 
discrete variables in the superconducting array. The or- 
thogonal superposition of V = ±1 states, ^(|+) — | — )), 
corresponds to the half-vortex inside the hole and has a 
much larger energy in the superconductive array. 

The degenerate ground states in the insulating array 
can be manipulated in the same way as in the super- 
conductor. Since physically these states correspond to 
the absence or presence of the half-vortex inside the hole 
the adiabatic change of local magnetic field that drags 
one half vortex across the system, flips the state of the 
system, providing us with the implementation of the op- 
erator t x acting on the state of the qubit. Analogously, 
motion of elementary charge 2e around the hole changes 
the relative phase of the states with and without half- 
vortex providing us with the operator r z . 

The signature of the topological insulator is the per- 
sistence of the trapped half flux inside the central hole 
(see Fig. 1) which can be observed by cycling magnetic 
field so as to drive the system back and forth between 



insulating and superconducting states. This trapping is 
especially striking in the insulator. Experimentally, this 
can be revealed by driving slowly the array into a super- 
conducting state and then measuring the phase difference 
between opposite points such as A and B in Fig. 1. In the 
state with a half vortex the phase difference is tt/2 + im 
while it is nn in the other state. The irn contribution is 
due to the usual vortices that get trapped in a big hole. 
This slow transformation can be achieved by changing 
the strength of weak links using the external magnetic 
field as a control parameter. The precise nature of the 
superconductive state is not essential because phase dif- 
ference 7r between points A and B can be interpreted as 
due to a full vortex trapped in a hole in a converntional 
superconductor or due to a 7r periodicity in a topologi- 
cal one which makes no essential difference. These flux 
trapping experiments are similar to the ones proposed 
for high-T c cuprates 8,16 with a number of important dif- 
ferences: the trapped flux is half of $o, the cycling does 
not involve temperature (avoiding problems with excita- 
tions) and the final state can be either conventional or 
topological superconductor. 
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